Wrinkling of microcapsules in shear flow 
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Elastic capsules can exhibit short wavelength wrinkling in external shear flow. We analyse this 
instability of the capsule shape and use the length scale separation between the capsule radius 
and the wrinkling wavelength to derive analytical results both for the threshold value of the shear 
rate and for the critical wave-length of the wrinkling. These results can be used to deduce elastic 
parameters from experiments. 
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The pronounced shape transformations of vesicles and 
microcapsules induced by hydrodynamic flow paradig- 
matically illustrate a main theme in microfluidics. The 
dynamical balance between fluid-induced viscous and 
mcmbranc-detcrmincd clastic stresses depends crucially 
and distinctively on the specific soft object immersed into 
the flow. In shear flow, fluid bilayer vesicles assume a sta- 
tionary tank-treading shape if there is no viscosity con- 
trast between the interior and exterior fluid 0. If the 
interior fluid or the membrane gets more viscous a tran- 
sition to a tumbling state can occur 0, 0, S S IE 0- 
While comprehensive experiments are still lacking, tank- 
treading has been observed both for vesicles in infinite 
shear flow and, in particular, for vesicles interacting 
with a rigid wall ^3 where a dynamical lift occurs 
[ill IT^ llil ll^ . For a comprehensive review of the dy- 
namical behaviour of soft capsules in shear flow, we refer 
the reader to the first two chapters of [T^ . 

In contrast to fluid vesicles, genuine microcapsules ex- 
hibit a finite shear rigidity since their membrane is phys- 
ically or chemically cross-linked. While shear elastic- 
ity prevents very large shape transformations such as a 
prolate-oblate transition in fluid vesicles , it also leads 
to qualitatively different instabilities in shear flow like the 
wrinkling instability first observed experimentally |l6l |. 
see Fig. ^ In fact, wrinkling is a ubiquitous phenomenon 
which can occur whenever a compressive force acts on a 
thin sheet either clamped at boundaries or attached to 
an elastic substrate as happens for the wrinkling of skin 

Eaili. 

The quantitative study of wrinkling is challenging since 
the usually large deformations of thin sheets lead to non- 
linear partial differential equations for the force balance. 
Still, rather general statements about wrinkling phenom- 
ena become possible, since the nonlinearity is mostly due 
to geometric effects and thus universal. Indeed, it has 
been pointed out by Cerda et al. in a series of intrigu- 
ing papers [13, lla, lis! that geometry puts such a strong 
constraint on the underlying physics, that wrinkles which 
form in very different contexts all obey the same scaling 
laws. They consequently identified as the necessary in- 
gredients for all wrinkling phenomena a thin sheet with a 
bending stiffness, an effective elastic foundation, and an 




FIG. 1: Wrinkling of a polysiloxane capsule (~ 343/im) in 
shear flow as indicated by the arrows (adapted from [T^])- 



imposed compressive strain [1^. These ingredients give 
rise to a rather general effective Hamiltonian description, 
where long wavelength wrinkles are punished energeti- 
cally due to the increasing longitudinal stretching strain, 
while short wavelength wrinkles increase the bending en- 
ergy. This leads to the selection of the optimal wrinkling 
wavelength. 

Since Cerda et al. considered universal scaling be- 
haviour of all wrinkling phenomena, conversely no con- 
clusions concerning the physical origin underlying the 
effective constants can be drawn from the scaling be- 
haviour alone. Considering the elastic capsule system 
in particular, identifying the effective elastic stiffness of 
the rather complex coupled fluid-membrane system is not 
trivial. Moreover, the Hamiltonian description is a priori 
no longer valid when the compressive forces acting on the 
membrane are nonconservative in nature, e.g. when aris- 
ing from a hydrodynamic flow. Here, dynamical effects 
have to be taken into account. In this letter we there- 
fore complement the approach discussed above by taking 
advantage of the short wavelength of the observed wrin- 
kles. We can thus identify a new mechanism by which 
short wavelength wrinkles are selected in a curved geom- 
etry. This approach can also be generalised to the full 
dynamic problem, although we have so far only taken 
into account the short-time dynamics. Moreover, the ex- 
act asymptotic limit allows us to make quantitative (i.e. 
non scaling argument) predictions on the onset and di- 
rection of wrinkles. 



2 



Assuming linear elasticity (and presupposing the cap- 
sule was initially stress-free) the work required to deform 
the membrane is given by the expression 

n[R] ^^JdA ^ATr(e)2 + 2^ie: e + K (^2H - ;| j ^ ■ 

. . 

Here A is the two-dimensional isotropic compressibility, 
fi the Lame coefficient corresponding to the shear stress, 
and K the bending modulus of the elastic material, for 
which we assume for simplicity a spontaneous curvature 
given by the initial capsule radius. The resulting forces 
can be calculated as functional derivatives of the energy 
with respect to the deformation f°' = — (5?i[R]/(5R. We 
now immerse this capsule of radius R into a fluid of vis- 
cosity 77 subject to the unperturbed shear flow 

V- = 7[(y/2)e, + {x/2)ey], (2) 

where 7 denotes the shear rate. (Here x and y refer to 
global Cartesian coordinates fixed in a laboratory frame.) 
Considering more general linear shear flow is straightfor- 
ward, but would distract from the wrinkling problem. 
For the system at hand the resulting Reynolds num- 
ber of the flow is negligible, so that the velocity field 
V and pressure p are governed by the Stokes equation 
77AV = Vp together with the incompressible continuity 
equation V • v = 0. With the definition of the hydro- 
dynamic stress cr = —pi + rj [Vv -f (Vv)-^] the Stokes 
equation can also be written as V • cr = 0. The hy- 
drodynamic force acting on the membrane is given by 
f** = [a -n]™'. The Stokes equations are accompanied 
by the asymptotic boundary condition v(r) v"^ for 
|r| 00, the force balance f°' + f ^ = at the mem- 
brane, and the kinematic condition that the membrane 
is convected by the fluid, vjp^^g rt>t)~ dtR{9, (f), t). 

These equations fully determine the evolution of the 
capsule shape. Unfortunately they are highly nonlinear 
and non-local, so that a solution is only possible numeri- 
cally or in certain asymptotic limits. Let us now suppose 
that we have found a stationary solution, and look for the 
stability of small perturbations i5r of the membrane posi- 
tion. From the experimental evidence [l^ it is expected 
that beyond a critical shear rate 7c the stationary defor- 
mation becomes unstable with respect to the formation of 
short wavelength wrinkles. For a first qualitative discus- 
sion, we consider a given configuration with known initial 
strain tensor eo and impose oscillatory wrinkles perpen- 
dicular to the membrane with the large wavevector /ce^-, 
leading to 5y = Sr cos{kx)ez. In this one-dimensional 
picture, the strain is simply the relative change in length 
of a line element with respect to a reference element. An 
initially straight line (a cut through a flat membrane) 
always increases in length when folded. The relative 
length change is of the order 6I1/I ~ {6r)'^k'^. This ef- 
fect is of second second order in the amplitude and grows 



quadratically with the wavenumber. For a curved mem- 
brane, there is an additional first order length change 
SI2/I ^ dr cos{kx) / R, arising from the fact that the lo- 
cal curvature radius of the line is periodically displaced 
outward and inward. Both contributions are added to 
the initial strain. For the curved membrane, the relevant 
component of the strain tensor thus changes to 

60 + [Srfk^ + Srcos{kx)/R. (3) 

Even though the curvature induced part is oscillatory and 
vanishes in the mean, it nevertheless gives rise to a posi- 
tive contribution in the elastic energy, which is quadratic 
in the strain. For large wavenumbers, the local curvature, 
which is of the order k'^SR, gives rise to the most impor- 
tant contribution to the total energy. After averaging 
over the fast oscillations, one obtains for the elastic en- 
ergy (ignoring pre-factors of the order one and the shear 
elasticity for a moment) 

n''^nf + J dA^{Sr f + Xeo{6r fe + K{6r)^k\ (4) 

This expression has formally the same form as the ef- 
fective Hamiltonian in |l8j . so that we can identify the 
bending stiffness k, the stress Aeo, and the effective elastic 
foundation stiffness X/R^. Before we discuss the latter 
quantity in more detail, let us briefly recapitulate how 
these terms lead to short wavelength wrinkling for suf- 
ficiently large compressive stress: For an initially com- 
pressed membrane (i.e. eo < 0), the compressive stress 
is partially canceled by the second order contribution 
((5r)^fc^. This lowers the modulus of the resulting stress 
and therefore decreases energy. This mechanism is also 
the driving force for Euler buckling of planar membranes 
and rods. For a curved membrane, there is an additional 
contribution that always serves to increase the elastic en- 
ergy. Only with a sufficiently large compressive pre-strain 
(eo < 0) one is able to overcome this energy cost. Since 
the energy gain is quadratic in the wavenumber, short 
wavelength wrinkles will overcome the cost associated 
with the curvature effect sooner for increasing | — eo|. 
However, very large wavenumbers are penalised by the 
bending energy. A simple discussion of the functional 
form Q thus reveals that gaining elastic energy by wrin- 
kling is possible provided that the initial compressive 
stress exceeds a threshold tq = A|eo| > ^ ^/Xk/R. 
Even though this scaling prediction is already implicit in 
[isl l the 1/R term shows the crucial influence of a curved 
geometry here. 

This critical stress must be provided by the imposed 
hydrodynamic flow. From a simple dimensional anal- 
ysis one can deduce that this stress is of the order 
T ~ TjjR. For the critical stress we find a critical shear 
rate 7c ^ ^/kA/ {r]R^) and a critical wrinkling wavenum- 
ber kc ~ {X/{kR^)Y/^ ^ l/(/ii?)^/^, where the second re- 
lation arises when we treat the membrane as an isotropic 



3 



shell of thickness h and bulk modulus E, i.e. A 
and Kr ^ Eh^. 



Eh 



This discussion reveals the physical origin of the effec- 
tive stiffness X/R^: Any normal deformation on a curved 
membrane leads to an additional strain contribution that 
increases the elastic energy to second order in the ampli- 
tude. This picture becomes more subtle when wc allow 
the membrane to relax tangentially: In plane displace- 
ments will reduce the second order strain to minimise 
its energy contribution. The two dimensional freedom in 
general does not allow to reduce the full stress tensor, 
which consists of three independent components, to zero. 
This becomes possible, however, for particular values of 
the material constants or the membrane geometry. If the 
shear modulus of the membrane vanishes, tangential de- 
formations can completely eliminate the isotropic strain. 
In such a membrane the effective stiffness consequently 
vanishes. This explains why wrinkling is not observed 
for fluid membranes. When the curvature perpendicu- 
lar to the wrinkle wavevector vanishes, it becomes again 
possible to eliminate the total additional stress. Thus 
the wrinkling phenomenon on a capsule necessarily needs 
both a curved geometry and finite values for the com- 
pressibility and shear modulus. 

The quantitative treatment is considerably more com- 
plicated. We not only have to express the change in the 
elastic energy due to wrinkling not necessarily perpen- 
dicular to the membrane on a quantitative level, but we 
also have to take the change in the hydrodynamic forces 
due to the deformed capsule into account, since change 
in the forcing linear in the amplitude acts on the same 
level as the quadratic change of the elastic energy. How- 
ever, we can neglect all these complications if we assume 
that the number of wrinkles on the membrane, kcR, is 
large. Formally, we are looking for the singular pertur- 
bation solution of the dynamic equations in the limit 
e = (K/(Ai?2))i/4 ^ l/{kcR) 0. When we compare 
the elastic forces with the change in the hydrodynamic 
forces in detail, we we find that the latter can be ignored 
as being of the order e. 

We now turn to the quantitative predictions. The wrin- 
kling pattern is assumed to be of the more general form 
i5r = Stq exp(i0), where the local wavevector is now given 
as the gradient of the phase k = V0. The amplitude 
Stq and the wavevector k are allowed to vary slowly on 
the membrane. The collection of the most important 
contributions to the chan ge in elastic energy is formally 
obtained via WKB theory |20j . We need to consider non- 
normal deformations of the membrane, which allows the 
membrane to reduce the stress even further. We will pub- 
lish the details of the calculation elsewhere and only note 
that the energy change due to the formation of wrinkles 



on a sphere is 



■k + K\k[ 



5r' 



(5) 

Here t denotes the tangential part of the elastic stress 
tensor, which is effectively a symmetric 2x2 matrix deter- 
mined by the external forces. When this matrix has a suf- 
ficiently large negative eigenvalue, energy can be gained 
from the formation of wrinkles. The force driving the 
folding is the negative derivative of with respect to 
5r and in this order linear in the amplitude. The elas- 
tic forces in turn induce a flow in the surrounding fluid 
that convect the membrane points. In general, this ef- 
fect would be non-local, leading to a coupling between 
the wrinkling dynamics at different points of the cap- 
sule. However, at zero order in the asymptotic short 
wavelength limit the hydrodynamic effects are local and 
diagonal in the coordinate system adapted to the mem- 
brane. We are thus led to the evolution equation 



dt{5r) 
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4/i A + ^ 
8ry|A:| i^A + 2/i 



k-T-k + K|/fcm<5r = ^, 
(6^ 

which implies a time scale r for the short time dynamics 
of each wrinkling mode. We take the time constant as 
a measure of the strength of the instability and assume 
that locally the fastest growing mode is dynamically se- 
lected. We can thus find the resulting folding pattern by 
maximising the r.h.s. of equation Note that the only 
input needed for this theory is the initial stress distribu- 
tion T, which is caused by the hydrodynamic flow. For 
general shapes, these data can come e.g. from numerical 
simulations such as 0. For the present purpose, we will 
use analytic results for the expansion of the stress to first 
order in the shear rate 7 [2l| . The stress distribution on a 
stationary sphere in linear shear flow (O is independent 
of the material constitutive law and reads in spherical 
coordinates, see Fig. |2 



T00 Tg^ 



.5 / sin 20 cos 20 cos 
^^2'' I cos 20 cos — sin20cos^6' 



(7) 

The wavevector field of the wrinkles is given by the direc- 
tion of the negative eigenvalues of (j?)). Even though these 
can be calculated analytically, they are given by rather 
involved expressions. However, one can see that the most 
negative stress eigenvalue — 57?/i?/2 appears along the 
circumference = 37r/4, 77r/4, which is at 7r/4 angle 
with the shear direction. The corresponding eigenvec- 
tor points perpendicular to the this circumference along 
the coordinate lines. Along this circumference the lo- 
cal compression is maximal, from which we can calculate 
the critical shear rate 7c and the corresponding critical 
wavenumber kc analytically as 
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(8) 



4 




30 




500 


25 




450 




400 


20 \ 




350 




300 






250 






200 




150 






100 






50 











FIG. 2: Wrinkling instability region for 7 = 11.5s~^ (left) and 
7 = 17.0s~^ (right), respectively. The elastic moduli of the 
membrane are X = fi = O.lNm"^, k = 1 ■ 10~^^Nm. In the 
black colored region the membrane is still stable with respect 
to wrinkle formations. Brighter colors indicate a larger growth 
rate of the wrinkles. Initially wrinkling occurs only in a small 
region around the circumference corresponding to the plane of 
maximal compressive stress. At higher shear rates this region 
grows, as does the growth rate of each wrinkling mode. The 
direction of the wrinkles along the lines of maximal tension 
are indicated schematically. 



These two expressions are our main quantitative results. 
The scaling is of course the same as predicted in [lij . 
However, we also predict quantitative prefactors that 
become exact in the large wavenumber limit. Above 
the critical shear rate jc wrinkles start to form in a fi- 
nite region starting uniformly along the circumference 
(j) = 37r/4, 77r/4. With increasing shear rates this wrin- 
kling region grows symmetrically towards the stagnation 
points of the flow located a,t 9 = tt/2, = 7r/4, 57r/4. 
The wrinkling direction is determined by the direction 
of the eigenvectors of the stress tensor. Along the initial 
wrinkling circumference, the folds form perpendicular to 
this line, i.e. in the direction of maximal tension, in ac- 
cordance with experimental findings (16j . 

Inser ting experimental value for the planar elastic 
moduli jig A ~ ^ ~ 0.1 Nm~^ and counting the number 
of wrinkles across the membrane (c.f. Figure^ we arrive 
at an estimate for the bending rigidity k ~ 1 • 10~^^ Nm, 
in accordance with a similar estimate made in [l^. Of 
course, due to the fourth power taken, the estimation 
for K (c.f. equation ((SJ) depends sensitively on the mea- 
sured wavenumber. For these values, together with the 
capsule radius R ~ 343 /xm and the viscosity of water, 
77 10"'^ Pa • s, the critical shear rate is predicted to be 
7c ~ lls~^. This value is in good agreement with the 
experimentally reported onset at jc — 4s~-^ For 
shear rates larger than 7c, the minimization of © is 
best performed numerically. The growth of the wrin- 
kling region for the numerical shear rates 7 = 11.5s~^ 
and 7 = 17.0s~^ is shown in Figure|21 

In conclusion, wc have presented an analytical theory 
for the wrinkling instability of a microcapsule in shear 
flow. This was accomplished by a systematic singular 
perturbation expansion in the bending rigidity. The only 
input needed, apart from the material constants, is the 



initial stress distribution. For quasispherical capsules all 
quantities can be calculated explicitly, leading to ana- 
lytical expressions of the shear rate and the wrinkling 
wavelength at the onset. This provides a method to de- 
duce the bending rigidity from shear experiments. Given 
the relatively large error in the experimental values of 
the compressibility the agreement of the predicted crit- 
ical shear rate (deduced from the observed wavelength) 
with experimental values is surprisingly good. 

Finally we remark that an initial stress on the mem- 
brane (e.g. due to osmotic pressure) affects the critical 
shear rate. While the effect of membrane prestress on the 
deformation of capsules has been discussed the ef- 
fect on the formation of wrinkles has not yet been consid- 
ered. Since wrinkling can only occur for membranes with 
shear stress resistance, fluid membrane capsules do not 
form wrinkles. Thus a study of viscoelastic membranes, 
which exhibit fluid or elastic behaviour depending on the 
timescale, seems worthwhile. 

Financial support of the DFG within the priority pro- 
gram "Nano- and Microfluidic" is acknowledged. 
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